Spin ordered phase transitions in 
isospin asymmetric nuclear matter 
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Spin polarized states in nuclear matter with Skyrme effective forces are studied on the base of a 
Fermi liquid theory for a wide range of isospin asymmetries and densities. There are a few possible 
scenarios of spin ordered phase transitions: (a) nuclear matter with SLy4 interaction undergoes at 
some critical density a phase transition to a spin polarized state with the oppositely directed spins of 
neutrons and protons; (b) for SkI5 interaction, a spin polarized state with the like-directed neutron 
and proton spins is formed; (c) nuclear matter with SkI3 interaction under increasing density, at 
first, undergoes a phase transition to the state with the opposite directions of neutron and proton 
spins, which goes over at larger density to the state with the same direction of nucleon spins. 
Spin polarized states at strong isospin asymmetry are accompanied by the long tails in the density 
profiles of the neutron spin polarization parameter near the critical density, if the energy gain of 
the transition from the nonpolarized state to a polarized one is the decreasing function of isospin 
asymmetry (SLy4 force). If the energy gain is increased with isospin asymmetry, there are no long 
tails in the density distribution of the neutron spin polarization parameter (SkI3, SkI5 forces). 

PACS numbers: 21.65.+I; 26.60.+C; 75.25.+z; 71.10.Ay 



The issue of spontaneous appearance of spin polarized 
states in nuclear matter is a topic of a great current in- 
terest due to its relevance in astrophysics. In particular, 
the scenarios of supernova explosion and cooling of neu- 
tron stars are essentially different, depending on whether 
nuclear matter is spin polarized or not. On the one hand, 
the models with the effective Skyrme and Gogny nucleon- 
nucleon (NN) interaction predict the occurrence of spin 
instability in nuclear matter at densities in the range from 
go to 6p n fo r different parametrizations of the NN poten- 
tial (go = 0.16 fm -3 is the nuclear saturation 
density). On the other hand, for the models with the 
realistic NN interaction, the ferromagnetic phase transi- 
tion seems to be suppressed up to densities well above 

go M-M- 

Here we continue the research of spin polarizabiiity of 
nuclear matter with the use of an effective NN interac- 
tion. The main objective is to study the possible sce- 
narios of spin ordered phase transitions in nuclear mat- 
ter with Skyrme forces, attracting parametrizations of 
a NN potential being relevant for calculations at strong 
isospin asymmetry and high density. In particular, we 
choose SLy4 effective interaction, constructed originally 
to reproduce the results of microscopic neutron matter 
calculations j23|. We utilize SkI3 and SkI5 parametriza- 
tions as well, giving a correct description of isotope shifts 
in neutron-rich medium and heavy nuclei p^ . As com- 
pared with the research of Ref. [12| with SLy4 and SLy5 
effective interactions, here we explore a wider domain of 
isospin asymmetries, including symmetric nuclear matter 
and neutron matter as limiting cases. Besides, we use 
also SkI3 and SkI5 parametrizations, that will allow us 
to study new scenarios of spin ordered phase transitions, 
not found in Ref. |T3 |. 

The basic formalism is presented in detail in Ref. ^2] • 
We are interested in studying spin polarized states with 



like-directed and oppositely directed spins of neutrons 
and protons. One should solve the self-consistent equa- 
tions for the coefficients £oo, £30, £o3j £33 m the expansion 
of the single particle energy in Pauli matrices in spin and 
isospin spaces 

Wp) = £ o(p) + £oo(p) - Moo, &o (p) = £30 (p) j (1) 
£03 (P) = £03(P) - M03, £33 (P) = £33 (p)- 

Here £o(p) is the free single particle spectrum, and 
£00 : £30: £03 : £33 are the Fermi liquid (FL) corrections to 
the free single particle spectrum, related to the normal 
FL amplitudes £7o(k), £73 (k) by formulas 



£~oo(p) = 2y^2u {k)f 00 (q), k 



p q 



(2) 



£30( P )-7^E C/ l( k )/3"( C l)' 

q 

£o 3 (p) = 7^E C/ 2(k)/o3(q), 



£33(p) = 7^E C/ 3(k)/33(q). 



The distribution functions /oo, /03, /30, /33, in turn, 
can be expressed in terms of the components £ of the 
single particle energy and satisfy the normalization con- 
ditions for the total density g n + g p = g, excess of 
neutrons over protons g n — g p = ag, ferromagnetic 
(FM) p-f — gi = Agjj and antiferromagnetic (AFM) 
(g n ij + g P i) — (g n [ + g p <\) = Agjj. spin order parameters, 
respectively (a being the isospin asymmetry parameter, 
Q\ = Qn\ + g P i and g± = g nl + g pl , with g nh g nl and 
ffpl ' 6pl being the neutron and proton number densities 
with spin up and spin down). The quantities of interest 
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are the neutron and proton spin polarization parameters 



n„ 



gn| ~ Qnj 
Qn 



n„ 



characterizing spin ordering in neutron and proton sub- 
systems. In the numerical solution of the self-consistent 
equations we utilize SLy4. SkI3 and SkI5 Skyrme forces. 
Note that the density dependence of the nuclear sym- 
metry energy, calculated with these Skyrme interactions, 
gives the neutron star models in a broad agreement with 
the observables j^. Another important constraint on 
the set of Skyrme force parameters can be obtained, if 
to consider expression for the effective mass of a neutron 
to* in totally spin polarized neutron matter 



to 



gm 



37 = 1 + ^2(l + ^ 2 ), 



(3) 



where too is the bare mass of a nucleoli. Eq. fol- 
lows from expressions for FL amplitudes in neutron mat- 
ter Since usually for Skyrme parametrizations ti < 
0, then we get the constraint xi < — 1, which guarantees 
the stability of totally polarized neutron matter at high 
densities ^3, 01 • The Skyrme parametrizations SLy4, 
SkI3 and SkI5 satisfy this condition. 

Fig. shows the density dependence of the neutron 
and proton spin polarization parameters at zero temper- 
ature for SLy4 force. The main qualitative feature is that 
for SLy4 force there are only solutions corresponding to 
the oppositely directed spins of neutrons and protons in 
a spin polarized state and there arc no solutions corre- 
sponding to the same direction of neutron and proton 
spins. The reason is that for SLy4 force the FL am- 
plitude Ui, determining spin-spin correlations, becomes 
more and more repulsive with the increase of nuclear mat- 
ter density, while the FL amplitude Us, describing spin- 
isospin correlations, becomes more and more attractive 
with density, leading to spin instability with the oppo- 
sitely directed spins of neutrons and protons at higher 
densities. The critical density of spin instability in sym- 
metric nuclear matter (a = 0), corresponding to AFM 
spin ordering (Agjj. ^ 0, Ag^ = 0), is g c ~ 0.33 fin -3 . 
It is less than the critical density of FM instability in 
neutron matter, g c w 0.59 fm~ 3 . Even small admixture 
of protons to neutron matter leads to the appearance of 
long tails in the density profiles of the neutron spin po- 
larization parameter near the transition point to a spin 
ordered state. As a consequence, a spin polarized state 
is formed much earlier in density than in pure neutron 
matter. 

As seen from Fig. ^d, for SkI5 force, oppositely to the 
case with SLy4 force, there are only solutions correspond- 
ing to the same direction of neutron and proton spins in a 
polarized state and there arc no solutions corresponding 
to their opposite directions. Comparing to the previous 
case, the FL amplitudes U\ and Us exchange their roles: 
the FL amplitude Us becomes more and more repulsive 
with the increase of nuclear matter density, while the FL 
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FIG. 1: (Color online) Neutron and proton spin polarization 
parameters as functions of density at zero temperature for (a) 
SLy4 force and (b) Ski 5 force. 



amplitude U\ becomes more and more attractive with 
density, leading to spin instability with the like-directed 
spins of neutrons and protons at higher densities. For 
SkI5 force, a phase transition to the FM spin state in neu- 
tron matter takes place at the critical density g c « 0.28 
fm~ 3 . It is less than the critical density of spin insta- 
bility in symmetric nuclear matter g c ss 0.43 fm -3 , cor- 
responding to FM spin ordering ( Ag>t| ^ 0, Ag^ = 0). 
An important peculiarity is that there are no long tails in 
the density profiles of the neutron spin polarization pa- 
rameter at large isospin asymmetry. Hence, in this case 
a small admixture of protons to neutron matter doesn't 
considerably change the critical density of spin instability 
and even leads to its increase. 

In Fig. |21 the total energies per nucleon of the spin 
ordered and nonpolarized states are compared at zero 
temperature for SLy4 and SkI5 forces. One can see that 
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FIG. 2: (Color online) Total energy per nucleon, measured 
from its value in the normal state, as a function of density at 
zero temperature for (a) SLy4 force and (b) SkI5 force. 



nuclear matter undergoes a phase transition to the state 
with the oppositely directed (SLy4 force) or like-directed 
(SkI5 force) spins of neutrons and protons at some criti- 
cal density, depending on isospin asymmetry. It is worth 
to note an important difference in isospin dependences 
for these two cases. For SLy4 interaction, the difference 
between the total energies of spin polarized and nonpo- 
larized states is largest at the given density for symmetric 
nuclear matter while for SkI5 interaction it is largest for 
neutron matter. This means that a phase transition in 
density to a spin polarized state will take place earlier 
in symmetric nuclear matter than in neutron matter for 
SLy4 force, while for SkI5 force, oppositely, it will occur 
earlier in neutron matter, as was mentioned above. 

Fig. shows the neutron and proton spin polariza- 
tion parameters as functions of density at zero tempera- 
ture for SkI3 force. There are two types of solutions of 
the self-consistent equations in symmetric nuclear mat- 
ter, corresponding to FM and AFM ordering of neu- 
tron and proton spins. Due to proximity of FL am- 



FIG. 3: (Color online) The dependences for SkI3 force: (a) 
Same as in Fig. 0with the legends of Fig. (b) same as in 
Fig. |21 Also the curves, corresponding to FM and AFM spin 
ordering in symmetric nuclear matter, are shown. 



plitudes Ui and U3, the respective critical densities are 
very close to each other (g c ps 0.910 fm -3 for FM order- 
ing and g c i=a 0.917 fm~ 3 for AFM ordering) and larger 
than the critical density of spin instability in neutron 
matter (g c w 0.37 fm~ 3 ). When some admixture of pro- 
tons is added to neutron matter, the last critical density 
is shifted to larger densities and a spin polarized state 
with the oppositely directed spins of neutrons and pro- 
tons appears. Under increasing density of nuclear matter, 
the neutron spin polarization continuously increases till 
all neutron spins will be aligned in the same direction. 
Protons, at first, become more polarized with density 
and their spin polarization is opposite to the spin polar- 
ization of neutrons. But, after reaching the maximum, 
spin polarization of protons decreases and at some criti- 
cal density spins of protons change direction, so that the 
spin ordered phase with the like-directed spins of neu- 
trons and protons is formed. Then, beyond the critical 
density, the spin polarization of protons is continuing to 
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increase until the totally polarized state with parallel or- 
dering of neutron and proton spins will be formed. Thus, 
for SkI3 force nuclear matter undergoes at some critical 
density a phase transition from the state with antiparal- 
lel ordering of neutron and proton spins to the state with 
parallel ordering of spins. With increasing isospin asym- 
metry, this critical density increases as well. Note that 
there are no long tails in the density profiles of the neu- 
tron spin polarization parameter at large asymmetries. 

In Fig. the total energies per nucleoli of spin or- 
dered and nonpolarized states are compared at zero tem- 
perature for SkI3 force. In symmetric nuclear matter, 
FM spin ordering is thermodynamically more preferable 
than AFM one. The energy gain of the transition from 
the nonpolarized state to a spin polarized state increases 
with isospin asymmetry at the given density, analogously 
to SkI5 force. Hence, a spin polarized state in neutron 
matter occurs earlier in density than in symmetric nu- 
clear matter, as was clarified above. Note that for SkI3 
force under increasing density, initially, the state with an- 
tiparallel ordering of neutron and proton spins appears 
without existence of long tails in the density profiles of 
the neutron spin polarization parameter at strong isospin 
asymmetry. This is in contrast with SLy4 force, for which 
the antiparallel ordering at a < 1 is characterized by 
the appearance of such long tails. Hence, the presence 
of long tails in the density profiles of the neutron spin 
polarization parameter doesn't relate to the antiparallel 
ordering of neutron and proton spins, but is associated 
with the decreasing dependence of the energy gain of the 
phase transition to a spin polarized state as a function 
of isospin asymmetry In this case, the critical density of 
spin instability in symmetric nuclear matter is less than 
the critical density in neutron matter. If these critical 
densities are substantially different, that is the case for 
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SLy4 force, then the long tails in the density profiles of 
the neutron spin polarization parameter appear. 

In summary, we have considered spin ordered phase 
transitions in nuclear matter with SLy4, SkI3 and SkI5 
Skyrme effective forces. It has been shown that asym- 
metric nuclear matter with SLy4 effective interaction un- 
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sition is characterized by the appearance of long tails in 
the density profiles of the neutron spin polarization pa- 
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metry. This means, that even small admixture of protons 
to neutron matter leads to the considerable shift of the 
critical density of spin instability to lower densities. The 
presence of such long tails is associated with the decreas- 
ing dependence of the energy gain of transition to a spin 
polarized state as a function of isospin asymmetry 

In nuclear matter with SkI5 effective interaction a spin 
polarized state with the like-directed spins of neutrons 
and protons is formed. For SkI3 force, nuclear matter un- 
der increasing density, at first, undergoes a phase transi- 
tion to a spin polarized state with the oppositely directed 
spins of neutrons and protons. Under further increasing 
density, spins of protons change their direction at some 
critical point and a phase transition from the state with 
antiparallel ordering to the state with parallel ordering of 
neutron and proton spins occurs. Since the energy gain 
of the transition from the nonpolarized state to a polar- 
ized one is increased with isospin asymmetry, there are no 
long tails in the density distribution of the neutron spin 
polarization parameter for SkI3 and SkI5 forces. Note 
that the obtained results may be of importance for the 
adequate description of thermal and magnetic properties 
of neutron stars. 
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